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It  is  shown  that  Taylor's  beampattern  for  a  continuous  aperture  can  be 
computed  analytically  without  Fourier  transforming  the  weighting  function 
itself,  thereby  achieving  economies  in  computational  effort  in  some  modeling 
situations.  A  short  Fortran  program  is  given.  An  approximate  formula  for  the 
half-power  beamwidth  is  derived.  It  is  pointed  out  that  the  Taylor  weighting 
function  can  be  negative  for  large  "n,  a  fact  that  does  not  seem  to  be  well 
known.  In  addition,  modification  of  Taylor's  design  to  force  the  weighting 
function  to  go  to  zero  as  a  power  o  of  distance  from  the  aperture  endpoints  is 
discussed.  For  a  =  1  and  a  =  2  this  results  in  an  increase  of  5%  and  10%, 
respectively,  in  the  beamwidth. 
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I.  INTRODUCTION 

This  Memorandum  is  a  review  of  Taylor's  original  weighting  function  for 
continuous  apertures.  It  is  presented  in  some  detail  in  Sections  II  and  III. 
It  is  shown  that  Taylor's  beampattern  and  weighting  function  can  be  computed 
easily  by  analytically  exact  formulas.  Taylor's  beampattern  turns  out  to  be 
the  product  of  a  rational  function  and  the  beampattern  of  a  uniformly  weighted 
aperture. 

Also  reviewed  is  a  modification  due  to  Rhodes  of  Taylor's  pattern  for  the 
purpose  of  forcing  the  weighting  function  to  go  zero  as  a  power  a  of  distance 
from  the  aperture  endpoints.  This  results  in  a  5 t  increase  in  beamwidth  over 
the  beamwidth  of  Taylor's  original  pattern  if  a  =  1,  and  a  10%  increase  if 
«  =  2  (for  TT  =  10;  see  below).  These  modifications  are  discussed  in  Section 
IV. 


Taylor's  original  paper  [1]  derives  a  symmetric  weighting  function  for  a 
continuous  aperture.  He  does  not  discuss  or  even  mention  its  use  for  arrays 
of  point  sensors.  His  method  is  essentially  an  ad  hoc,  but  intuitively 
sensible,  procedure  which  blends  together  the  desirable  characteristics  of 
uniform  weighting  and  the  van  der  Maas  weighting  into  one  weighting  design. 

The  blending  is  accomplished  by  careful  specification  of  the  beampattern 
nulls.  The  various  si  delobe  levels  do  not  enter  the  method's  derivation.  In 
other  words,  the  si  delobes  are  whatever  they  turn  out  to  be  after 
specification  of  the  nulls. 

It  is  often  said  that  Taylor  weighting  makes  the  first  few  sidelobes  near 
the  mainlobe  nearly  flat;  that  is,  all  "near-in"  sidelobes  have  essentially 
the  same  amplitude.  Thjs  statement  is  erroneous.  See  Figure  1,  for  example, 
where  the  9  sidelobes  (n  =  10)  nearest  the  mainlobe  would  all  be  at  -20  dB  if 
the  statement  were  true.  Instead,  the  first  sidelobe  is  at  -20  dB  and  the 
ninth  sidelobe  is  at  (roughly)  -25  dB. 

It  is  a  useful  fact  that  the  beampattern  corresponding  to  Taylor 
weighting  can  be  computed  analytically,  without  Fourier  transforming  the 
weighting  function.  This  can  be  seen  from  Taylor's  original  discussion  [1], 
which  is  reviewed  in  this  Memorandum.  Taylor's  original  notation  is  retained 
here.  Appendix  A  gives  a  FORTRAN  program  which  computes  the  beanpattern 
and/or  the  weighting  function  using  the  analytical  formulas  developed  below. 

In  addition,  it  computes  the  exact  half-power  beamwidth. 


The  aperture  is  assumed  to  lie  on  the  p-interval  from  -it  to  +ir.  The 
weighting  function  g(p)  is  related  to  the  far-field  beampattern  F(z)  by 


g(p)  e1zp  dp. 


(1) 


Taylor  assumes  throughout  that  g(p)  is  a  real  even  function.  Consequently, 
F(z)  is  also  an  even  function  of  z. 


It  is  well  known  that 
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where  eg  =  1  and  em  =  2  for  m  >  1.  In  other  words,  knowledge  of  the 
integer  samples  of  F(z)  implies  knowledge  of  F(z)  everywhere.  A  very 
different  representation  of  F(z)  is  the  infinite  product 


where  {zi,  Z2,  ...}  is  a  complete  list  of  all  the  positive  zeros  of  F(z). 

It  is  an  interesting  mathematical  fact  that  these  zeros  must  all  lie  on  the 
real  z-axis.  For  example,  uniform  weighting  g(p)  =  l/(2w)  gives 


F(z) 


sin  ifz  , 
irZ 


whose  positive  nulls  are  {1,2,3,  From  (3),  then. 


sin  ttz 
irz 


.77 


n  =  1 


(1 - j)» 


(4) 


a  well  known  identity  dating  back  at  least  to  Euler's  time  (circa  1750), 


By  means  of  his  choice  of  nulls  { zn  }  in  the  representation  (3)  of  F(z), 
Taylor  sought  a  beampattern  which  had  a  flat  envelope  near  the  mainbeam  and, 
for  large  z,  an  asymptotic  6  dB/octave  decay  rate.  He  also  sought  by  this 
same  means  a  physically  realizable  aperture  to  approximate  the  physically 
unrealizable  ideal  van  der  Maas  function.  (It  is  unrealizable  because  of  the 
presence  of  delta  function  spikes  at  the  aperture  end-points,p  =  ±  ir.)  Taylor 
found  a  set  of  nulls  which  came  close  to  attaining  his  first  objective  and 
which  did  attain  his  second  objective.  The  next  section  is  a  description  of 
Taylor’s  nulls. 


II.  TAYLOR'S  NULL  SPECIFICATION 


Taylor  specifies  the  nulls  z«  of  his  beampattern,  starting  with  n  =  n, 
to  be  exactly  the  same  as  those  of  the  uniform  weighting  function;  that  is, 


zn  =  n  for  n  =  n,  n  +  1,  ...  (5) 

The  positive  integer  n’  is  a  free  parameter  which  can  be  chosen  as  desired. 
Note  that  n  =  1  gives  exactly  uniform  shading.  Note  also  that  the  null  list 
(5)  guarantees  a  6  dB/octave  asymptotic  decay  rate  as  z-»  «.  .  (This  follows 
from  (8)  below.) 
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To  complete  the  list  of  positive  nulls  for  his  beampattern,  Taylor 
selects  (when  n  >  1)  the  "near-in"  nulls  to  be 


z_  =  CT 


iG 


(n  - 


for  n  =  1,2, 


n  -  1, 


(6) 


where 


a  =  n  /  ▼  A2  +  (n  -  y) 

A  --i  In  (R  *  Vr2  -  1) 

R  »  h',{/U 

S  =  maximum  sidelobe  level  (in  dB). 

This  choice  for  the  first  ”n— 1  nulls  may  seem  mysterious  at  first  glance,  but 
it  is  a  choice  based  on  the  ideal  van  der  Maas  function  [2],  defined  by 

Fq(z,  A)  =  cos  ir^z2  -  A2,  A>0. 

It  is  an  interesting  mathematical  fact  that  among  all  beampattern  functions 
F(z)  such  that 

(a)  F(z)  has  a  Fourier  transform  vanishing  outside  the  aperture  -it  to  +ir 

(b)  |  F ( z) | <_  1  for  \z  |  >  A, 

the  one  with  the  maximum  possible  value  at  z=0  is  the  van  der  Maas  function 
F(z)  =  Fq(z,A).  The  positive  nulls  of  Fq(z,A)  are 


=  W  +  (n  -  j)2,  n  -  1,  2,  3,  ... 

Comparison  of  these  nulls  with  Taylor's  ad  hoc  null  specification  (fc)  shows 
that  Taylor's  nulls  are  related  to  the  van  der  Maas  nulls  by  a  dilation  factor 
a.  The  factor  o  is  chosen  to  be  slightly  larger  than  unity  to  compensate  for 
the  6  dB/octave  decay  of  the  beampattern  for  z  >  n.  Note  that  n  =  «>  gives 
exactly  the  van  der  Maas  beampattern. 

III.  TAYLOR'S  BEAMPATTERN  AND  WEIGHTING  FUNCTION. 

Taylor's  beampattern  can  now  be  expressed,  using  (3),  as 
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The  last  expression  in  (7)  can  be  rewritten,  using  (4),  to  give 


sin  irz 

irZ 


(8) 


In  this  expression,  limits  must  be  taken  whenever  z  =  0,  1,  2,  3,  ...,  n-1  to 
avoid  the  indeterminate  form  0/0.  See  Appendix  B.  Note  that  Taylor's 
beampattern  is  identically  the  product  of  a  rational  function  (of  degree “n— 1 
in  z2)  and  the  beampattern  of  the  uniformly  weighted  aperture,  sin(nz)/irz. 

It  is  clear  that  Taylor's  beampattern  can  be  computed  analytically  from 
(8)  without  computing  the  weighting  function  at  all.  The  representation  (2) 
of  F(z)  is 


n-1 

F(z)  2ir  y''  «m 
m=0 

since  F(n)  =  0  for  n  >  n.  This  is  not  as  efficient  as  using  (8).  However,  it 
does  yield  an  efficient  way  to  compute  the  weighting  function  g(p).  By 
Fourier  transforming  it  term  by  term  and  using  the  fact  that  F(m)  =  F(-m).we 
get 


F  (m) 


(sin  ir 

~ iTz 


(z-m) 
z-m) 


sin  ir(z+m) 
n(z+m) 


•) 


(9) 


g(p)  = 


n-1 

F (0 )  +  2^J F(m)  cos  mp 
n=l 


.  |P|  <  w 


(10) 


This  is  the  (spatial)  Fourier  series  of  Taylor's  weighting  function.  By 
computing  once  and  for  all  the  constants  F(0),  F(l),  ...,  F(n-l)  using  (8), 
the  series  (9)  can  be  an  efficient  formula  for  computation. 


The  beamwidth  measured  between  the  first  nulls  is  (from  (6)  with  n=l) 


BWNULL 


2 


where  a  and  A  are  given  as  above.  An  exact  formula  for  the  half-power 
beamwidth  is  not  available.  Table  2  gives  half-power  beamwidths  that  were 
computed  numerically  (using  a  general  purpose  subroutine  in  [3,  Chapter  7]  ). 
More  useful  perhaps  is  the  following  approximate  formula  for  the  half-power 
beamwidth 
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BW 


3dB 


n-1 


2 

it 

t r 


♦*E 


n-l  \a2  (A2+  (n -h  ) 


-1/2 


To  prove  (11),  note  that  the  asymptotic  expansion 


F(z)  -  1  - 


n-1 


a2  (A2  +  (n-  -i)2)  ^^n^ 

n=l  n=n 


z  ,  z->  0 


follows  immediately  from  (7).  Since 


(ID 


n-1 


n-1 


•2  V  1 

■  r-Z-7* 

n=l 


we  have 


n-1 


F(z)  =  1  - 


*E 


~2 —  I o  • 


tfV2(A2  +  ("-7))  " 

Setting  F(z)  =  1/2  and  solving  for  z  gives  (11). 

The  accuracy  of  (11)  is  good  in  two  limiting  cases.  As  rv-x*',  o-»  1 
and  (11)  becomes 

/  °° 


BW 


3dB 


^E 


t?  a2  +  '"-T> 

The  exact  answer  for  the  van  der  Maas  function  is 


-1/2 


(12) 


BW3dB  "  2 


A2  -  ^  arc  cosh  (-|  cosh  uA)^ 


1 1/2 


and  a  comparison  with  (12)  is  given  in  the  last  row  in  Table  3.  Similarly, 
for  rT  =  1 ,  the  sum  in  (11)  vanishes  and 


.  2  VF 

=  -  =  i«  . 


BWojr  =  *"■  =  1.103  radians 

oats  ir 
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which  is  within  10  percent  of  the  correct  answer  of  BI^b  =  1.207  radians 
for  the  uniformly  weighted  aperture. 

Table  3  gives  the  relative  error  between  the  approximation  (11)  and  the 
exact  half-power  beamwidth  for  the  same  entries  as  in  Table  2.  It  may  be 
concluded  from  Table  3  that 

(a)  the  approximation  (11)  is  always  on  the  low  side  of  the  exact 
half-power  beamwidth,  and 

(b)  the  correction  required  to  make  (11)  exact  is  a  constant  factor 
which  depends  strongly  on  the  specified  sidelobe  level  and  very 
weakly  on  TT. 

Consequently,  a  suitable  correction  factor  depending  only  on  specified 
sidelobe  level  would  make  (11)  very  accurate. 

The  Taylor  weighting  function  need  not  always  be  a  positive  function. 
The  best  way  to  show  this  is  by  example.  Consider  the  case  n  =  100  and  a 
sidelobe  level  of  S  =  -20  dB.  The  weighting  function  is  slightly  negative 
just  inside  the  aperture  endpoints  (for  p  =±3.078761,  for  example,  Taylor's 
weight  is  -.005519929).  See  Figure  2.  The  Taylor  function  in  practice  is 
nearly  always  positive  for  smaller  values  of  n. 


n _ 

-lOdB  (A  =  .578) 

-20dB  (A=  .953) 

-30dB  ( A=  1.32) 

-40dB  (A=l. 69) 

5 

1.0475 

1.3264 

1.5526 

1.7323 

10 

1.0009 

1.2818 

1.5220 

1.7262 

15 

.9851 

1.2641 

1.5051 

1.7126 

20 

.9  771 

1.2548 

1.4954 

1.7036 

25 

.9724 

1.2491 

1.4892 

1.6975 

30 

.9692 

1.2452 

1.4849 

1.6932 

100 

.9581 

1.2313 

1.4691 

1.6761 

00 

.9533 

1.2252 

1.4619 

1.6680 

Table 

2.  Exact  Taylor 

half-power  beamwidths 

FT _ 

-lOdB  (A=. 578 

-20dB  ( A=i9  53) 

-30dB  (A=l. 32) 

-40dB  (A=l. 69) 

5 

7.67  % 

9.98  % 

11.4% 

12.2% 

10 

7.57  % 

9.91  % 

11.3% 

12.2% 

15 

7.55  % 

9.90  % 

11.3°% 

12.2% 

20 

7.54% 

9.89  % 

11.3% 

12.2% 

25 

7.55% 

9.89  % 

11.3% 

12.2  % 

30 

7.54% 

9.89  % 

11.3% 

12.2% 

100 

7.55  % 

9.88  % 

11.3% 

12.1% 

00 

7.54  % 

9.88  % 

11.3% 

12.1% 

Table  3.  Relative  error  of  approximation  (11)  to 
Taylor  half-power  beamwidths. 
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IV.  MODIFICATIONS  OF  TAYLOR  WEIGHTING 

Rhodes  [4,5]  shows  that  the  Taylor  weighting  function  g(p)  can  be  made  to 
go  to  zero  as  any  power  a  >  -1  of  distance  from  the  aperture  endpoints  by 
altering  the  position  of  the  nulls  in  Taylor's  function  F(z).  The  general 
design  technique  can  be  viewed  as  an  extension  of  certain  ideas  in  Taylor's 
original  paper  [1],  using  mathematical  methods  developed  by  Rhodes.  The  most 
important  cases  are 

1.  a  =  0,  which  is  exactly  Taylor's  original  case;  F(z)  decays 

asymptotically  at  6  dB  per  octave. 

2.  a  =  1,  for  which  the  weighting  function  goes  to  zero  linearly  at  the 

aperture  endpoints;  F ( z)  decays  asymptotically  at  12  dB  per  octave. 

3.  a  =  2,  for  which  the  weighting  function  goes  to  zero  quadratically  at 

the  aperture  endpoints;  F(z)  decays  asymptotically  at  18  dB  per  octave. 

The  cases  a  =  1  and  a  =  2  are  given  explicitly  below,  after  giving  the  method 
for  any  «  >  -1. 

A  theoretically  significant  criticism  of  Rhodes'  work  is  that  he  does  not 
prove  that  his  technique  is  mathematically  correct.  The  available  theory  (due 
to  Paley  and  Wiener,  and  to  Levinson)  provides  a  proof  only  for  -1/2  <  «  <  1/2 
and  a=l.  As  Rhodes  states  [5],  "it  is  not  unreasonable  to  expect  that  the 
general  theory"  is  valid  for  all  a  >  -1.  In  any  event,  we  can  proceed  to 
develop  the  method  for  all  a  >  -1  in  a  purely  formal  way,  ignoring  a 
theoretical  question  which  may  in  the  end  not  be  of  any  practical  importance. 
Taylor's  original  method  is,  after  all,  an  ad  hoc  technique  and  so  is  Rhodes’ 
generalization  of  it. 

Rhodes'  development  retains  the  integer  n  as  the  breakpoint  between  the 
near-in  nulls,  which  are  dilated  versions  of  van  der  Maas'  nulls,  and  the 
outer  nulls,  which  force  the  asymptotic  decay  rate  for  F(z)  to  be  6(1+  o)  dB 
per  octave.  Consequently,  in  the  limit  as  TT-*00,  the  van  der  Maas  function  is 
again  obtained  for  all  a  >  -1,  just  as  in  Taylor's  original  design  a  =  0. 

This  means  that  the  desired  behavior  of  the  weighting  function  at  the  aperture 
endpoints  is  confined  to  small  neighborhoods  of  the  aperture  endpoints  for 
larger  n.  In  other  words,  the  weighting  function  changes  rapidly  just  inside 
the  aperture  endpoints  for  large  n. 

The  development  in  [4]  is  brief  and  only  the  case  a  =  1  is  given  in  any 
detail.  His  later  paper  [5]  gives  enough  detail  to  carry  out  the  general 
development  for  a  >  -1.  This  requires  the  identity,  valid  for  a  >  -1, 


(1  +  a/2) 


r(z  +  1  -  a/2)  sin  n(z  -  a/2) 
r(Z  +  1  +  a/2)  ir(z  -  a/2) 


(13) 


It  is  proved  as  follows.  A  special  case  (zi  =  z?  =  a/2  and  Z3  =  z)  of  a 
result  in  [6,  Equ.  1.3(4)]  gives 
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r (z  +  a/2 


urn 

j  r  (-2 


Tm 


1  - 


(n  +  £) 


Dividing  by  the  first  term  in  the  infinite  product,  and  then  using  the 
recurrence  formula  [6,  Equ.  6.1.15]  and  the  reflection  formula  [6,  Equ. 
6.1.17]  whenever  necessary,  gives 


T  (z)  =  _ (a/2)2  r2(a/2) _ 

a  [(a/2)2  -  Z2]  r (a/ 2  +  z)  r(a/2  -  z) 

_ r  2(1  +  a/2) 

r(l  +  z  +  a/2)  r (l  -  z  +  a/2) 
p2(-|  +  1)  p(z  -  -|)  1 

’r(i  +  f +  z)  p(z  -f)  rd  -  (z  -|)) 

r2(f  M)p(z-f)  sinir(z-f) 

r(z  +  f  +  i)  u  (14) 

Multiplying  and  dividing  by  z  -  (a/2)  on  the  right  hand  side  of  the  last 
equation  yields  (13).  (We  note  that  above  is  given  without  proof  by  Taylor 
[1,  Equ.  (29)].) 


Rhodes  defines  the  general  Taylor  pattern,  Fa(z),  for  a  >  -1  to  be 


(15) 


where 


a 


a 


(n  +f) 


(16) 


and  A  is  the  same  as  given  above  (just  after  (6)).  Note  that  for  a  *  0  the 
function  Fo(z)  is  exactly  Taylor's  original  function  F(z).  The  analog  of 
(8)  for  general  a  is 
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n-1 


'1  - 


-7TI 


c\{ A2  +  (n  -\)) 


n-1 


1  - 


(n  t-f)2 


(17) 


as  is  clear  from  (15)  and  (13).  As  ,  the  rational  function  of  degree 

n-1  in  z2  in  (17)  approaches  a  constant  and  Ta(z)  is  asymptotic  to  a 
constant  (depending  only  on  o)  times  l/|z|l+a.  The  asymptotic  decay  rate  of 
Fa(z)  is  therefore  6(1  +  a)  dB  per  octave.  In  addition,  the  asymptotic 
decay  rate  means  that  Fa(z)  has  a  Fourier  transform  vanishing  outside  the 
aperture  [-ir,  ir]  for  every  a  >  -1  and  n  >  1. 

For  a  >  -1  define  the  "sampling  functions" 


(z) 


r(z  -  |  +  1)  sin  n(z  -  -|) 
T(z  +f)  *(z2  -  (n  +  f)2) 


(18) 


where 


I(-l)n(2n  +  a) (n  +  a)/n!  ,  if  a  4  0 
1  ,  if  a  =  0,  n  =  0 

( —1 ) n2  ,  if  o=0,  n=  1,2,3,... 

Each  g(®)(z)  is  an  even  function  of  z.  These  functions  are  essentially  the 

Lagrange  interpolating  functions  for  the  points  ^(n  +  a/2);  n  =  0,1,2,...}. 
More  precisely,  the  only  nulls  of  G*  (z)  are  of  the  form  ±(n  +  al 2)  and, 
furthermore. 


G^a) (±(m 


if 

if 


m  =  n 
m  4  n. 


(19) 


The  functions  G^a)(z)  are  derived  using  methods  due  originally  to  Paley 
and  Wiener. 


The  open  theoretical  question  mentioned  earlier  in  this  section  concerns 
the  completeness  of  the  sampling  functions  (18)  with  respect  to  all, even 
aperture-limited  functions.  As  stated  already,  it  is  known  that  G^®/  are 

complete  for  -1/2  <  a  <  1/2  and  for  <*  =  1.  For, other  values  of  a  >  -1,  nothing 
is  known.  Proceeding  on  the  assumption  that  G^®/  are  complete  for  all 

a  >  -1,  it  follows  Fa(z)  can  be. expanded  in  the  form 

oo 

F„<z> 

n=0 
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for  some  constants  an.  From  (19)  it  follows  that  an  =  Fa(n  +  (a/2))  for 
all  n.  From  (15)  it  follows  that  an  =  0  for  n  >  n.  Therefore, 


n  -  1 


F„<2> 


-Z)  Gn“)(2) 


(20) 


n  =  0 


which  generalizes  (9). 


Denote  by  ga(p)  the  weighting  function  corresponding  to  Fa(z).  As 
just  stated,  ga(p)  vanishes  outside  the  aperture  [-ir,  it].  Question:  Does 
9«(p)  9°  zero  as  the  Power  a  >  -1  of  distance  from  the  aperture 
endpoints  ?  Taylor  [1]  proves  that  any  even  function  with  this  endpoint 
behavior  has  a  Fourier  transform  whose  nulls  are  asymptotic  to  *(n  +  (a/2)), 
but  he  does  NOT  prove  the  converse.  Consequently,  although  ga(p)  is  even 
and  has  a  Fourier  transform  with  the  proper  null  locations,  this  is  not 
necessarily  sufficient  to  answer  the  question  in  the  affirmative.  However, 
taking  the  term  by  term  Fourier  transform  of  (20)  gives  the  expansion 


n  -  1 


9(.)<P)  *  S  F«(n  +  f>  Hfl°)(p) 


n  =  0 
■ 

where  for  n  *  0,  1,2,  ... 


(21) 


H<a)(p) 


I  ■ 

_  a  V  ^  (a) 

(2  cos  f)  2^  (~D 


r  -  0 


TrwTT • 


0 


|P|  <  * 


|P|>  T 


(22) 


where  eo  =  1  and  er  =  2  for  r  >  0,  and 

II  for  all  a,  if  k  =  0 

(a  +  1)  ...  (a  +  k  -  1)  for  all  a  ,  if  k  >  0. 

Since  each  of  the  functions  H^“)(p)  has  the  correct  endpoint  behavior, 
ga(p)  must  also  have  this  same  behavior. 


Just  as  in  Taylor's  original  case,  both  the  aperture  function  gtt(p) 
and  the  beampattern  function  Fa(z)  can  be  computed  independently  of  each 
other  using  the  analytically  exact  formulas  (21)  and  (17),  respectively,  for 
any  a  >  -1.  Appropriate  approximations  near  the  points  ±(n  +  a/2)  analogous 
to  those  developed  in  Appendix  B  for  a  =  0,  are  necessary  for  computing  Fa(z) 
using  (17).  Developing  these  approximations  should  not  present  any 
mathematical  difficulties. 

The  three  cases  a  =  0,  1,  2  are  now  given  explicitly.  Fortran  programs 
implementing  these  three  cases  should  be  easy  to  write.  The  sampling 
functions  are 
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G(°>(Z)  .  (.1)"  in.  ^  *  |i" 

it  (z  -  n  ) 

(23) 

G(1>(Z)  -  (-l)n+1  (2"+l)  cos  « 

n  *  (z2  -  (n+-l)2) 

(24) 

g(2)(z)  =  (-i)n+1  2<nn^  s1n  ,Z  -  . 

n  ir  z  (z2  -  (n-1)2) 

(25) 

the  corresponding  aperture  basis  functions  are 

Hp°^(p)  =  t-  cos  np 

(26) 

H^(p)  =  cos  (n+-|)p 

(27) 

H^2^(p)  =  (-l)n  +  cos  (n+l)p  . 

(28) 

It  should  be  noted  that  (23)  and  (26)  are,  within  a  scale  factor,  identical  to 
(9)  and  (10),  respectively.  In  all  cases,  the  aperture  function  ga(p)  is 
computed  from  (21).  Consequently,  the  only  potential  difficulty  is  computing 
the  constants  Fa(n  +  a/2)  for  is  n  =  0,  1,  ...»  n  -1.  Fortunately,  for 
o=0,  1,  and  2,  these  constants  are  easy  to  compute  using  (17)  since  the 
following  identities  hold: 


Tx(z) 

T2(z) 


sin  irz 


irz 

cos  irZ 


sin  irz 

irz  (1-Z2)  ’ 


(29) 

(30) 

(31) 


The  price  paid  for  the  desired  end  effects  is  an  increase  in  the 
beamwidth  over  the  beamwidth  of  Taylor's  original  weighting  function, 
beamwidth  measured  to  the  first  null  is,  for  all  a  >  -1, 


BW 


NULL 


+ 


1 


The 


where  aa  is  given  by  (16)  above.  This  gives  exactly,  for  fixed  n, 

%ULL  (for  an>  »>  ^  g 

BWNULL  ^for  °  °0  2n 

This  means  that  for  TT  =  10,  the  beamwidth  measured  between  nulls  is  5%  larger 
for  o=l  and  10%  larger  for  o  =  2  than  for  Taylor's  original  o  =  0 
beampattern.  It  is  anticipated  that  approximately  the  same  percentage 
increases  occur  in  the  half-power  beamwidths. 
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A  different  modification  to  Taylor's  nulls  can  be  utilized  to  produce 
asymmetric  beampatterns  using  complex  valued  aperture  functions  g(p).  This  is 
described  in  [8]  for  an  application  in  radar  to  minimize  ground  clutter.  The 
magnitude  of  the  aperture  function  turns  out  to  be  an  even  function,  while  the 
phase  of  the  aperture  function  turns  out  to  be  odd. 

V.  CONCLUSIONS 

Taylor  weighting  can  be  modified  to  force  the  weighting  function  to  go 
zero  as  any  power  o  >  -1  of  distance  from  the  aperture  endpoints.  Taylor's 
original  weighting  (a  =  0)  results  in  a  pedestal,  while  for  a  -  1  the 
weighting  function  goes  to  zero  linearly  as  in  a  cosine  window^and  for  a  =  2 
the  weighting  function  behaves  like  a  cosine-squared  window  at  the  aperture 
endpoints.  The  endpoint  effect  is  achieved  for  a  modest  increase  in  the 
mainlobe  beamwidth. 
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Appendix.  /\ , 

COMPUTE  TAYLOR'S  CONTINUOUS  Sri  ADI  nG  FUNCTION  AND  TRANSFER  FUNCTION 
T  N  PUT  REQUIREMENT:  1  .LE.  NBAR 

:  THE  SPATIAL  APERTURE  r,lES  FROM  -PI  TO  +  PI 
ARGUMENT  DEFINITIONS: 

x  =  abscissas  for  saropiina  the  snadlm  function 
s  =  shading  function  values 

ns  =  number  of  s  samples;  none  computed  if  ns  =  0 
k  =  abscissas  for  semolina  tne  transfer  function 
t  =  transfer  function  values 

nk  =  number  of  f  samples;  none  commuted  if  nk  =  0 
nbar=  the  first  nbar-1  zeros  of  f  are  tnose  of  van  der  Maas 
dp  =  sidelobe  level  in  db  of  limit ina  nolph-Cheovshev  array 
fir,  =  coefficients  of  the  shading  function  cosine  series 
b,v3ab  =  -3  dB  beamwidth;  not  computed  if  DW3do  is  set  to  -1. 
DIMENSION  LIMITS: 

x  and  s  must  be  dimensioned  at  least  max ( ns , +1 ) 
k  and  f  must  be  dimensioned  at  least  max ( nk ,  +  1 ) 
fm  must  be  dimensioned  at  least  near 
TEChNICAL  NOTES: 

NBAR=1  GIVES  THE  UNIFORM  SHADING  FUNCTION 
N e A R  =  1 N F I W I T Y  GIVES  PQLPH-CKEB YSH£ V  SHADING 
FIPST  BEAMPATTERN  HULL  =  SIGMA  *  SORT C a ** 2 +  .  7 5 ) 

THE  COSINE  SERIES  FOR  S  HAS  DEGREE  EXACTLY  NBAR-1 
PROGRAMMER:  R.  L.  5TPEIT,  NUSC,  DECEMBER  21,  1984. 

LAST  REVISION:  JANUARY  11,  1985 

subroutine  taylor(x,s,ns,k,t,nk,nbar,db,fm,pw3db) 

double  precision  x(l),sCl),k(l),f(l),db,fm(l),xpt,a , sigma , a, Pi, q  , 
zeroin,ap,bp,tol,ow3dp 
data  pi/3. 1415926535897980/ 
a=10.0d0**abs(db/20.0ri0) 
a=dlogCa+sqrtta*a-l,Odh))/pi 
slama=nbar/sqrt(a*a+(nbar-,5d0)*(nbar-.5d0)) 
nbar 1 =nbar-l 
fmC  1  )  =  1  .  oQ0 
i ± ( nbar . eq . 1 ) go  to  15 
do  10  i=2, nbar 
x  p  t  =  i  -  1 

f m ( i ) =g ( xpt , a , nbarl , sigma) 
continue 

if (ns. le.9)go  to  25 
do  20  1=1 , ns 
s (i )=0.0d0 

if(aDS(x(i)).gt.(pi+l.d-7))go  to  20 

s ( l ) =g ( x ( i ) , fm , nbar) 

continue 

i r  ( nk . le . 0 ) go  to  35 
do  30  1=1, nk 

f(i)=q(k(i)  ,a,noarl ,siama) 
continue 

it(bw3db.lt.0.0d0)go  to  40 
ap=0.0d0 

bp=sigma*sart(a*a+.25dO) 

tol=1.0d-10 

bw3db=2.0d0*zeroin(ap,bp,q,toi,a,nbarl,sigma) 

continue 

return 

end 

double  precision  function  g ( z , a , nbar 1 , s ioma ) 
double  precision  pi , Diz , zn , z , a , sigma 
data  pi /3. 14l59265358979d0/ 
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do  10  k  =  o , near  1 

l£Cabs(z-k).lt.l.0d-4)oo  to  3  0 
10  continue 

a=l . OaO 

i£(nbarl .ea.Ojao  to  21 
do20n=i,noarl 

zn=(.z/siama)**2/(a*a+(n-.5d0)**2) 
q=a* ( 1 .  OdO-zn )/ ( l.0d0-(z/n)**2) 

20  continue 

21  piz=oi*z 

c.=  (sln(piz)/piz)*q 
return 

30  itCk.gt.OJgo  to  50 

o=l . OaO 

1 f ( nbar 1 . eg . 0 ) ao  to  41 
do  40  n  =  l  ,noar 1 

zn=tz/sinma )**?/( a*a+(n-.5d0)**2) 
a  =  u*  ( 1 . Od  o-zn  j / ( 1 ,0a0-(z/n)**2) 

40  continue 

41  niz=pi*z 

a=(1.0d0-plz*piz*(1.0dn-piz*oiz/20.d0)/6.0dQ)*q 

return 

50  a=l . OaO 

do  bO  n= 1 , nba r 1 
1  £ ( n . ea . k ) ao  to  60 
zn=(z/sigrna)**2/(a*a+(n-,5d0)**2) 
o  =  g*  ( 1 .OdO-zn j/(l .0d0-(z/nl**2) 

60  continue 

zn=(z/sigina)**2/(a*a  +  (k-.bdC)**2) 
q=a*(l .OdO-zn) 
niz=pi*  (z-k ) 

a=(1.0d0-piz*plz*(1.0dO-piz*piz/20.d0)/6.0d0)*q 

o  =  a*(-l)**(l<+l)*(K/(z+z*z/K )) 

return 

end 

double  precision  function  a(p,fm,nbar) 
double  precision  p , twooi , f m ( 1 ) 
data  twopi/6. 293  165307  l  705>36477dO/ 
a=0 . OaO 

i f ( nba r . eq .  1 ) go  to  20 
do  in  i=2,nbar 
q=a+fm(i)*cos(  C i- 1 ) *p ) 

10  continue 

20  a=CfmCl)+g+g)/twopi 

return 
end 
c 

c  Compute  a  zero  of  a  real  function  f  in  the  interval  [ax,  nxj , 

c  Double  precision  version  of  program  on  po .  164-166  of  "Computer 

c  Methods  for  Mathematical  Computations,"  by  G.E..  Forsythe, 

c  v.A.  Malcolm,  and  C.B.  Moler,  Prentice-Hall,  1077,  put  slightly 

c  altered  for  use  in  computing  Taylor's  half-power  oeamwidth. 

double  precision  function  zer oin (ax , bx , £ , tnl , adum , near  1 , siqma ) 
double  precision  ax , ox , f , tol , adum , sigma 
do u ole  precision  a,b,c,d,e,ens,fa,fb,fc,toll,xm,p,a,r,s 
eps=l . ndO 

10  eps=eps/2,0d0 

to 11=1 .OdO+eps 
if (toll .ot.l.OdOJgo  to  10 
a  =  ax 
b  =  ox 
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fa=f(a,adum,nbarl,siama)-.5d0 
fo  =  r(o,adum,nDart  ,  s lama)-. 5d0 
20  c  =  a 

f  c  =  f  a 
d  =  o-a 
e  =  d 

30  1 1 Cabs ( f c J  .ge.aos (f o) )qo  to  40 

a  =  b 
b  =  c 
c=a 
f  d  =  £b 
f  p  =  r  c 
f.  c  =  £a 

40  toli=2.0dG*eps*ab5(p)+.5d0*toi 

xm= . 5d0  * ( c-o ) 

i£(abs(xm)  .le.tolllgn  to  90 
lf(tb.eq,0.0d0)go  to  90 
1  £  Cabs (e ). It .toll Jqo  to  70 
i f ( abs ( fa ) . le . abs ( f o ) ) do  to  70 
if (a.ne.c)qo  to  50 
s=f b/f a 
o  =  2 . OdO  *xm*s 
d=l .OaO-s 
oo  to  60 
50  q=f a/f c 

r=f b/f c 
s=f b/£a 

o=s*(2.0d0*xm*o*(q-r)-(b-a)*(r-1.0d0)) 

q=Cq-l ,0d0)*cr-1.0d0)*(s-1.0rt0) 

60  1 f ( P.gt . o . OdO ) q=-q 

p=aDS(o) 

i£(  C2.0d0*p)  .  ge .  (3.0d0*xo5lIq-abs(toll*q)  1  lao  to  70 
if (p.ge.abs(0.5d0*e*q) )go  to  70 
e=d 
d  =  n/q 
go  to  80 

70  d  =  x  m 

e=d 

80  a  =  b 

f  a  =  f  b 

i  f  ( aos  (  d )  .qt,tolUb=b+d 
ifCabs(d).le.toll) b=b+s ign ( toi 1 ,xm) 
fb=£(b,adum,nbarl,siama)-.5d0 
i r ( ( f b* C £c/abs Cf c) ) ) . gt . 0 . OdO) ao  to  20 
go  to  30 

90  zeroln=D 

return 
end 
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Appendix  8.  Calculation  of  F(z) 
For  some  small  number  e  >  0,  say  e  =  10-4,  define 


s(z,  k) 


sin  irz  i 

TtTkT  J  I 

(-l)k  £l--^(z-k)2  + 


z-k | >e>0 

■rar  (z_k)4 


Now,  if  | z— k|  >  e  for  k=  0,1,2,...,  n-1,  compute  F(z)  exactly  as  in  (8). 
fz— k|  <_  e  for  k=0,  then  compute 


If 


If  |z— k  |<_  e  for  k=l,  ...»  n-1,  then  compute 


Use  of  these  formulae  eliminates  all  indeterminate  0/0  forms  that  arise  during 
actual  computation  using  (8). 


■ 
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